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Abstract 

We present a search for the decay B + — > t + v t in 88.9 x 10 6 T(4S) decays recorded with the 
BABAR detector at the SLAC i?-Factory. A sample of semi-exclusive, semi-leptonic B decays 
{B~ — * D°£~i>£X), where X is either a photon, 7r°, or nothing, is used to identify the daughter 
particles that are associated with the other B meson in each event. These particles are searched 
for evidence of a one-prong leptonic B + — > t + v t decay. For this data sample we set a preliminary 
upper limit on the branching fraction of B(B + — > t + u t ) < 4.9 x 10 -4 at the 90% confidence level. 
This result is then combined with a statistically independent BABAR search for B + — > r + i/ T to give 
a combined preliminary limit of B{B + — > t + z^ t ) < 4.1 x 10 . 
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1 Introduction 



In the Standard Model, the purely leptonic decay B + — > r + z/ T 1 proceeds via quark annihilation 
into a W boson (Fig. Its amplitude is thus proportional to the product of the -B-decay constant, 
fs and the quark-mixing-matrix element V ub - The branching fraction is given by: 

ol2 



Gi 



T V) 



,m B m 7 
~8^ 



1 



m: 



T B+fB I Ke6 



(1) 



where we have set h = c = 1, Gp is the Fermi constant, V^, is the quark mixing matrix element 
HI E! /b is the i? + meson decay constant, describing the overlap of the quark wave- functions 
inside the meson, t b + is the B + lifetime and m B and m T are the B + meson and r masses. This 
expression is entirely analogous to that for pion decay. Physics beyond the Standard Model, such 
as SUSY, could enhance this process, as through the introduction of a charged Higgs 

Current theoretical values for fg have large uncertainty, and purely leptonic decays of the B + 
meson may be the only clean experimental method of measuring fg precisely. Given measurements 
of Vub from semi- leptonic processes such as B — > iriv, fs could be extracted from the measurement 
of the B + — > t + v t branching fraction. In addition, by combining the branching fraction measure- 
ment with results from B mixing the ratio |Kift|/|^td| can be extracted from B(B + — > t + v t )/ Am, 
where Am is the mass difference between the heavy and light neutral B meson states. 




Figure 1: Purely leptonic B decay B + — > £ + v^ proceeds via quark annihilation into a W boson. 

The decay amplitude is proportional to the lepton mass and decay to the lighter leptons is 
severely suppressed. This mode is thus the most promising for discovery at existing experiments. 
However, challenges such as the large missing momentum from several neutrinos make the signature 
for B + —s- t + v t less distinctive than for other leptonic modes. 

The Standard Model estimate of this branching fraction is 7.5 x 10 -5 , given the PDG 2002 
values for f B , \V ub \, and the B + lifetime: f B = 198MeV, V ub = 0.0036, t b + = 1.674 ps HJ. 

Purely leptonic B decays have not yet been observed. CLEO and experiments at LEP [£j [Zj |H1 
have searched for this process and set limits on the branching fraction at the 90% confidence level 
(CL). The most stringent limit on B(B + — > t + v t ) comes from the L3 experiment: 

B{B + t + v t ) < 5.7 x 10" 4 at the 90% CL (2) 

2 The BABAR Detector and Dataset 



The data used in this analysis were collected with the BABAR detector at the PEP-II storage 
ring. The sample corresponds to an integrated luminosity of 81.9 fb _1 at the T(4S) resonance 

1 Charge-conjugate modes are implied throughout this paper. The signal B will always be denoted as a B + decay 
while the semi-leptonic B will be denoted as a B~ to avoid confusion. 
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(on-resonance) and 9.58 fb~ taken 40MeV below the T(4S) resonance (off-resonance). The on- 
resonance sample consists of about 88.9 million BB pairs. The collider is operated with asymmetric 
beam energies, producing a boost of flj ~ 0.56 of the T(45 l ) along the collision axis. 

The BABAR detector is described elsewhere jOJ. Charged-particle momentum, direction, and 
energy loss (^§) measurements are performed by a five-layer double-sided silicon vertex tracker 
(SVT) and a 40-layer drift chamber (DCH) which are immersed in the field of a 1.5-T super- 
conducting solenoid. Charged particle identification is performed by a detector of internally- 
reflected Cherenkov light (DIRC). The energies of electrons and photons are measured by the 
electromagnetic calorimeter (EMC) consisting of an array of 6580 CsI(Tl) crystals. The flux return 
is instrumented with resistive plate chambers (IFR) to detect the passage of muons and neutral 
hadrons. 

A GEANT4-based Monte Carlo (MC) simulation is used to model the signal efficiency 
and the physics backgrounds. Simulation samples equivalent to approximately three times the 
accumulated data were used to model BB events, and samples equivalent to approximately 1.5 
times the accumulated data were used to model e + e~ — > uu, dd, ss, cc, and t + t~ events. 

3 Analysis Method 
3.1 Event Selection 

The decay B + — > r + z/ T and subsequent one-prong leptonic decays r + — ► e + v e V T and r + — > fi + u fJ T' T 
lead to the production of a single, visible charged particle and missing energy carried by neutrinos. 
Therefore, any remaining neutral energy or charged tracks originate from the other B. The search 
for B + — > t + u t proceeds by reconstructing one of the B mesons in each event in the semi- leptonic 
topology B —¥ D°£~Dj)X (henceforth referred to as the semi-leptonic side or semi-leptonic B), 
where X could be a photon, neutral pion, or nothing. The remainder of the event is searched for 
evidence of the one-prong leptonic r + decay modes. 

The choice of a semi-exclusive, semi-leptonic B decay means that the B meson cannot be fully 
reconstructed, since there is a neutrino missing and excited neutral D states are not reconstructed, 
potentially leaving unassigned neutral energy in the event. The inherent cleanliness of the events 
(small multiplicity, little unassigned neutral energy), however, is a powerful constraint that compen- 
sates for the choice of reconstruction mode. This method has been used in a search for B + — > K + vv 

m. 

Events are pre-selected by requiring that their net charge is zero, that R2 < 0.9 (where R2 
is the normalized second Fox- Wolfram moment J2]), and that they have large missing mass 
(M m i ss > 1.0 GeV). We select our semi-leptonic side by first reconstructing D° candidates in 
one of four modes: D° -> K~ir + , D° -> R-tt+tt , D° -> K-tt+tt+tt-, and D° -> if°7r + 7r". The 
Kg is reconstructed only in the charged mode K® — > 7r + 7r~. The D meson is required to have a 
momentum in the center-of-mass (T(4S)) frame p* D0 > 0.5 GeV, and its daughters are required to 
meet particle identification criteria consistent with the particle hypothesis. 

We pair the D° candidates with leptons whose momentum in the center-of-mass frame meets 
the requirement p\ > 1.3 GeV. The D°-lepton pairs are required to converge at a common vertex; 
these pairs are referred to as D£ candidates. If there are multiple D£ candidates in an event we 
select the candidate whose reconstructed D mass is closest to the average reconstructed mass for 
its decay mode. If the best candidate can be paired with a pion to make a D* + l~ candidate, then 
the event is consistent with being a B°B° event and is rejected. 
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Figure 2: The D° mass distribution after making event preselections and requiring that events have 
at least one B~ — ► D°£~P(X candidate. Background MC combines BB (excluding B + — > r + z/ T 
decays) and continuum MC and is normalized to the data area, while the normalization of the 
signal simulation is arbitrary. Due to the large size of this sample, the data error bars are too small 
to see in the figure. 

The mass distribution for these D° candidates is shown in Fig. El There is a difference between 
the mean of data and the mean of MC because the MC generator uses a D° mass that differs 
from the PDG value by 1 MeV. The width of distribution is typically narrower in MC by 0.2 MeV 
(K~tt + ) to 4 MeV (i^ _ 7r + 7r°) due to energy and momentum resolution differences in data and MC. 
To account for differences in mean position and peak width, we require that the reconstructed D° 
mass be within 3a of the fitted mean for the mode under consideration. The widths of the peaks 
vary by D° mode, from 3.8 MeV for D° K~tt + tt~'k + to 10.4 MeV for D° -► R-tt+tt . 

A restriction is also placed on cos 9b,d£i where 9g di is computed from B and Dl momenta and 
energy in the center-of-mass frame under the assumption that only a neutrino is missing, 

fl _ (2E B E m -m B -m 2 m ) 
COS 9 B ,Dt - ||^ — i • (a) 

2\PB\\PDe\ 

Since the semideptonic B cannot be completely reconstructed its energy, E B , and momentum, \pb\, 
are derived from the beam energies. The variable cos 9b, D£ can take on values outside (—1, 1) for 
events in which particles besides the neutrino have been missed in the semideptonic decay. We 
allow for the feed-down of excited neutral D states into the analysis by applying an asymmetric 
cut of —2.5 < cos 9b, De < 1-1. This variable is shown in Fig. |21 Differences in the semideptonic 
side between data and simulation are later corrected (Sec. 14. 2 j) in the overall efficiency estimate. 

Once the semideptonic side has been reconstructed, its daughter particles are excluded from 
the event. Any neutral or charged object that is not associated with the semideptonic side is 
referred to as belonging to the signal side. The event is then required to have only one remaining 
charged particle with a small impact parameter in the (r, <f>) plane. This track is required to fail 
kaon identification and be identified either as a muon or an electron. In addition, r + r _ continuum 
events are rejected by making requirements on the angle of the track with respect to the event 
thrust axis (| cos#-^| < 0.9) and the minimum mass capable of being made from any three tracks 
in the event {M™ in > 1.5 GeV). In general, continuum events are sharply peaked toward 1.0 in 
| cos#y^|; additionally, r + r~ events tend to have a value of M™ in peaked below the r mass. 
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Figure 3: The quantity cos Ob d£, computed from the energies and momenta of the D£ and the B 
meson assuming only a neutrino is missing. Background MC combines BB (excluding B + — > t + v t 
decays) and continuum MC and is normalized to the data area, while the normalization of the 
signal simulation is arbitrary. The arrows represent the boundaries of the selection region between 
—2.5 and 1.1. Due to the large size of this sample, the data error bars are too small to see in the 
figure. 

The variable E\ e ^, the neutral energy remaining in the calorimeter after excluding neutrals 
from the semi-leptonic B, provides good discrimination between signal and background (Fig. 0J. 
For true signal events, this quantity is peaked toward zero energy while for background it peaks 
near 2GeV. We make a loose restriction of -Ei e ft < 1.0 GeV in preparation for the final step of the 
analysis, which is to fit the -EWt distribution and extract the signal and background contributions. 

To minimize experimenter bias, the data were not studied in the signal-populated region of -Eieft 
until systematic checks were performed using control samples in the region of E\ e f t (E\ e f t > 0.5 GeV) 
which does not contain a significant signal yield. The region 0.5 < .Eieft < 1-0 GeV is defined to be 
the neutral energy sideband, while the region E\ e f t < 0.5 GeV is referred to as the signal region. 

3.2 Selection Efficiency 

After applying all selection criteria we obtain a signal efficiency in the Monte Carlo simulation of 
(6.99±0.31) x 10~ 4 , where the error is purely from the statistics of the signal MC. This number does 
not include corrections for differences between data and simulation, which will be described in Sec. 
|1J The number of background events expected for £i e ft < 1.0 GeV, determined from simulation, is 
269.1 ± 9.6, where the error is purely statistical. 

3.3 The Limit-Setting Procedure 

We fit -EWt using an extended maximum likelihood function; from this fit we extract the contri- 
butions to the distribution from signal and background events. The likelihood function is of the 
following form: 

n 

c = <,-(*«.+/*) -Q Fs (E lem ) + fi b F b (E Mtl )) (4) 

i=l 
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Figure 4: ^i e ft, the neutral energy remaining in the calorimeter after excluding neutrals associated 
with the semi-leptonic side. In the above distribution all analysis selection criteria are applied. The 
data are not shown below 0.5 GeV. The normalization of the signal MC is arbitrary. 

where fj, s is the number of signal events, /ib is the number of background events (both of which are 
obtained by fitting the selected data sample) and n is the total number of events in the selected 
sample. The functions F s (E\ e ft) and Fb(Ei e &) are the Probability Density Functions (PDFs) which 
are obtained from the E\ e [ t distributions in signal and background MC. The nominal models for 
signal and background are shown in Fig. |SJ 

The E\ c f t distribution in the signal MC has several salient features. The peak at 200 MeV arises 
from the unassigned photon or 7r° daughter of a D*° on the semi-leptonic side. The peak at zero 
energy arises from events without remaining neutral energy, and above zero the spectrum falls 
exponentially under the peak from the D*° neutrals, the result of neutrals in the MC arising from 
beam background and hadronic split-offs. 

The background PDF is a second-order polynomial rising monotonically between GeV and 
1.0 GeV. An estimate of the background from data is determined by extrapolating the data £i e ft 
sideband into the signal region using the background PDF. The result is a background estimate of 
273.6 ± 19.2 events. 

To determine the 90% confidence limit, denoted for the measured value of [i s we use a 
modified frequentist approach based on that used by the LEP experiments in the Higgs search |14j . 
This approach constructs an estimator based on the fitted signal yield, /j,{ ltted . In all following dis- 
cussions, toy MC refers to E\ e f t spectra obtained by randomly sampling the signal and background 
E\ e ft PDFs and treating the samples as a data set. The limit-setting method is outlined as follows: 

1. Toy MC experiments are generated with our expected number of background events and a 
signal hypothesis, //g npu * = ^expected_ j^g^ ^cy experiment is fitted for the signal yield, [i' s . 
The collection of fi' s forms a distribution. We denote the normalized fitted signal distribution 
for the original signal hypothesis as n™ put ). 

2. Toy MC is also generated with a signal hypothesis of zero events and with the expected number 
of background events. Again, each data sample is fitted and we obtain the distribution of the 
fitted signal yield. We denote the normalized distribution of the fitted signal yield in the null 
signal hypothesis as Db(n' s ). 

3. For a specific data sample, a fit is performed using the likelihood function to obtain the number 
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Figure 5: The signal PDF (left) fitted to E\ e f t from the signal MC and the background PDF (right) 
fitted to -Eieft from the background MC. All selection criteria are applied to the events in signal 
and background MC. The normalization of the signal MC is arbitrary and the normalization of the 
background MC is fixed to the integrated luminosity. 



of signal events in the sample, /j,{ ltted . The confidence level (CL) for this measurement, for a 
signal hypothesis jj,^ x P ected ^ [ s given by: 

^fitted £ d I 

nT i fitted expected^ _ I-oo ^ s+b{^' s ] /J,^ xpec e )dfl s 

4. The 90% confidence limit corresponding to this fitted signal yield, fi{ ttted , is defined by the 
requirement: 

CL(jjt{ med ,lif ) ) = 1-0.90 (6) 

A scan is performed over signal hypotheses until the minimum signal hypothesis which satisfies 
the requirement is determined. This signal hypothesis is The relationship between 
and [i{ ttted is given in Fig. H3 

We use the background-only hypothesis to determine the nominal sensitivity for the analysis. 
The nominal background PDF is used to generate 10,000 toy MC experiments, each with a total 
number of events Poisson-distributed around the nominal background expectation of 269 ± 9.6 
events. The distribution so generated for each experiment is fitted with the likelihood function, 
and the curve in Fig. [H] is used to obtain the corresponding (j,® . The distribution of fig for these 
toy experiments is shown in Fig. El and the median of this distribution provides our nominal 
sensitivity, fx". = 15.4. 

A test of the method is performed using two independent samples of GEANT4 MC, each with a 
total luminosity equal to that of the data set. Each simulated experiment contains only background 
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Figure 6: The relationship between fj,® and fj,{ ltted is shown in the left figure. In the right figure 
is the distribution of fj^° obtained from fitting 10,000 toy MC experiments generated using the 
null signal hypothesis. The black line indicates the median of this distribution, which is taken as 
the nominal sensitivity of the analysis = 15.4). The red and blue lines indicate the results of 
applying the analysis to two independent, background-only GEANT4 MC experiments. 

events and has all analysis selection criteria applied. The results obtained from fitting these two 
experiments, fj,® = 13.3 and [i® = 19.4 are in good agreement with the nominal sensitivity and 
are illustrated in Fig. |HI The results of these GEANT4 MC experiments are used in Sec. 14.11 to 
explore potential systematic effects on the limit-setting procedure. 

The strength of this method is that it takes full advantage of the information contained in the 
shape of the signal and background and is insensitive to the absolute background normalization. 
It relies on the accuracy of the simulation, in particular the simulation of the calorimeter for low 
neutral cluster energy. The systematic effects on the limit-setting procedure, as well as on the 
signal efficiency, are discussed in the next section. 

4 Evaluation of Systematic Uncertainties 

In this section we discuss the systematic uncertainties affecting the limit-setting procedure and the 
calculation of the efficiency. 

4.1 Systematic Effects for Limit-Setting 

TableHlists the sources of systematic uncertainty and their effects on the overall expected sensitivity 
of a pair of independent "experiments" derived from the background MC sample. Since MC is used 
to generate the PDFs for signal and background, effects of concern are those which might cause the 
true background or signal distributions to differ from the simulation. For each source of potential 
systematic uncertainty listed below, a CL curve is calculated (as in Fig. EJ) by generating toy MC 
with a systematically modified model and fitting with the nominal model. 

To explore the potential differences between data and MC when validating the background 



14 



PDF, several control samples have been developed by performing the semi-leptonic reconstruction 
on data and MC and applying selection criteria to the remainder of the event: 

• (B~ — > D°£~r>gX) + (anything): The B~ — > D G i~vgX reconstruction algorithm is applied 
and no restrictions are placed on the signal side. 

• (B° — > D* + £~v) + (1 track) : The semi-leptonic reconstruction algorithm is applied, but 
semi-leptonic candidates are kept where the best candidate is reconstructed as a D* + £~v. 
We then require only one track remaining in the event. This sample has a topology similar 
to the signal sample. 

• (B~ — > D®£~ v^X) + (2 tracks): The B~ — > D^l^v^X reconstruction algorithm is applied 
and the signal side is required to contain exactly two charged tracks. This requirement 
enforces a mis-reconstructed signal side, since total event charge is non-zero. 

The difference between data and simulation for each of these samples is used to construct a 
correction function for the PDFs. The PDFs are multiplied by the correction function and re- 
normalized. The upper limit is recalculated using the corrected PDFs to generate the toy MC 
samples discussed in Sec. while fitting is still performed with the nominal model. The results of 
this study are given in Tabled 

Motivated by the data/MC yield disagreement in the control samples, we explored the effects 
of a significant difference between the actual background and the background expectation. Given 
the expected background from MC (269.1 ±9.6 events), we varied the input background hypothesis 
to our toy MC by ±10 and ±20 events and used the nominal likelihood function to generate CL 
curves for each case. The result was a change in the nominal sensitivity of fi®° of no more than 0.5 
events (Table ^) . Therefore we concluded that if the actual background in data is different than 
the expectation from MC, the resulting change in the limit is small. 

The choice of background parametrization has also been explored as a possible effect on the 
limit-setting procedure. Different models (second-order polynomial, Gaussian tail, KEYS |13j . 
histogram) were used to represent the background E\ e f t spectrum and then generate toy MC. These 
toy samples were again fitted with the nominal model and the variation in the confidence limit 
curve was explored. The results are given in Tabled Limit curves calculated from this variation 
in the background model yielded a lower limit on the signal yield, suggesting that the second-order 
polynomial leads to a conservative limit. 

A feature at or near zero energy in the E\ e ft spectrum, unmodelled in the Monte Carlo, is 
another possible effect in the background that could fake a signal. Such a feature is suggested by 
studying the B° — > D* + £~v + 1 track control sample, which suggested a maximal difference near 
zero energy in data and MC of a factor of two. To explore the effect of such a feature, toy MC is 
generated with the background PDF enhanced by a factor of two near zero. This toy MC is then 
fit using the nominal signal and background models to study the variation in the confidence limit. 
The results are given in Tabled 

To validate the signal PDF, a control sample of B~ — > D°£~i>eX,B + — > D*°£ + u^ is used to 
check the difference between data and simulation after making all semi-leptonic and signal-side 
requirements. This sample has a semi-leptonic side as in the signal analysis, and a signal side 
which is also reconstructed as an independent B + — ► D°£ + i/£X. We then construct a subsample of 
these events where a photon or ir° could be associated with the D° on the signal side to make a 
D*°. Having attempted to account for the neutral from D*° decay on the signal side, the neutral 
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energy spectrum should resemble that for signal events. However, this subsample has more neutral 
energy present from hadronic splitoffs, and the correct photon or 7r° is not always associated with 
the D° to make the D . Therefore, a comparison is made between data and the corresponding MC 
sample and not between data and signal MC. The sample suggests that there may be a difference 
of ~ 20MeV in the overall position of the peak near 200 MeV. The effects of such a difference are 
given in Tabled 

Table 1: Sources of potential systematic uncertainty have been explored and their effects on the 
expectation for fj,® observed. The simulated experiments come from two independent samples 
of full MC, each with luminosity equivalent to the data luminosity. Despite the many potential 
differences between data and simulation, the variation in the upper limit is small copared to the 
variation of limit expected for the zero signal hypothesis. 



Systematic Variation 


Control Sample used 


90 


90 




to guide the variation 


Simulated 


Simulated 






Experiment 1 


Experiment 2 


None (Nominal Model) 


None 


19.4 


13.3 


Data/ Simulation 


B — > JJ u i VjiX + anything 


ion 

18. 


12.4 


Shape Difference 


B° _» D* + £~v + 1 track 


18.6 


13.0 




B~ -» D°e-i> e X + 2 tracks 


17.8 


12.2 


Incorrect Background Estimate 








-20 Events 




19.0 


12.0 


-10 Events 


All 


19.3 


13.1 


+10 Events 




19.7 


13.6 


+20 Events 




19.8 


13.7 


Gaussian Tail Background Model 


None 


17.1 


12.0 


KEYS Background Model 


None 


17.4 


12.0 


Histogram Background Model 


None 


17.2 


11.8 


Factor of 2 








enhancement in 


B° _» D* + £~v + 1 track 


16.5 


11.5 


background near 








£ left = 0.0GeV 








+20 MeV Shift in 








Signal PDF Gaussian 


B- -> D°£-v e X,B+ -» D*°l+v e 


20.2 


13.9 


-20 MeV Shift in 








Signal PDF Gaussian 


B- -> D°£-i? e X,B + -» D*°f + Vi 


18.9 


13.0 



In almost all cases that were investigated, the nominal model yields the most conservative CL 
curve. Additionally, the variation in the upper limit expectation is small compared to the overall 
variation of the limit for the zero signal hypothesis. We therefore use the nominal model to set the 
limit, having determined that the systematic effects mentioned above have negligible impact on the 
result. 
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Table 2: The corrections which are applied in this analysis on the signal efficiency, as well as the 
error associated with each correction. 



Source 


Correction Factor 


Correction Error 


B-Counting 


1.0 


0.011 


bemi-Leptonic 
Reconstruction Efficiency 


0.958 


0.051 


Neutral Energy 


1.0 


0.025 


Tracking Efficiency 
(signal track only) 


1.0 


0.008 


Particle ID 
(signal track only) 


0.836 


0.042 


Total 


0.801 


0.063 



4.2 Systematic Effects for Signal Efficiency 

The efficiency calculated in Sec. I3.2l needs to be corrected for differences in data and MC. The effects 
of several systematic uncertainty have been evaluated. The results are expressed in a correction 
factor and a systematic error associated with the correction. The most prominent sources of 
systematic uncertainties are listed in Table 

The semi-leptonic reconstruction efficiency correction is evaluated by using a sample of N 2 
double-reconstruction events (B~ — > D°£~i>iX, B + — > D°i + ui>X) and a sample of N\ single- 
reconstruction events in the same data set, where the single-reconstruction sample excludes double- 
reconstructed events. To isolate the effects of the semi-leptonic reconstruction, we remove the event 
preselection requirements which are unrelated to the selection of the semi-leptonic decay. The ratio 
of the sizes of these two samples is given by: 

Ni 2(1 -e) 

N 2 e { ' 

where e is the probability of selecting a second B meson in the mode B~ — > D°£~i>gX having 
already selected one B meson in that mode. Equation Q is solved for e, and the ratio of e in data 
and simulation is taken as a correction on the efficiency, 

£ -^- = 0.958 ±0.051. (8) 

£MC 

The correction of the signal efficiency due to particle identification (PID) is performed as follows. 
Using control samples of electrons and muons from data, the performance of the kaon veto and 
lepton identification is evaluated as a function of lab momentum (p) and polar angle (9). This 
performance translates into a weight for an identified track with a given momentum. Integrating 
this weight over the momentum spectrum of the signal MC for true one-prong leptonic r + decays 
provides the efficiency of the selection in data. 

The efficiency in the MC is evaluated by comparing the true identity of a t + daughter to the 
reconstructed identity. The ratio of the data efficiency to the MC efficiency is taken as the signal 
efficiency correction from particle identification. The statistical uncertainties on the signal MC 
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sample size and PID weights translate into a systematic error on the correction. This error is 
dominated by the signal MC sample size. 

We apply the total correction listed in Table |2] to the signal efficiency from MC to obtain the 
signal efficiency, (5.60 ± 0.25(stat.) ± 0A4(syst.)) x 1CP 4 . The uncertainty on each correction is 
propagated through to a systematic error on the signal efficiency, while the statistical uncertainty 
derives from the signal MC sample size. 

Systematic uncertainty is incorporated into the branching fraction limit as follows: for each 
branching fraction hypothesis, the corresponding signal hypothesis is determined by sampling 
10,000 times from a Gaussian distribution about the mean of the signal efficiency, where the width 
of the Gaussian is given by the systematic error on the signal efficiency. From each signal hypothe- 
sis a toy MC experiment is generated and fitted using the likelihood function. This is repeated for 
many branching fraction hypotheses and the CL curve, including systematic error, is calculated. 

5 Results 

We apply all selection criteria in the analysis and fit the -EWt spectrum in the data below 1.0 GeV 
using the nominal likelihood function. The result of the fit is given in Table El and shown in Fig. 
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Figure 7: The E\ e f t spectrum in data from zero to 1.0 GeV and the fit of the likelihood function to 
the spectrum. 



Table 3: The results of the fit of the likelihood function to data is given below. 



Fit Parameter 


Fitted Value 


, fitted 


10.9 ± 7.5 


fitted 


258 ± 17 



The background yield of 258 =L 17 is in agreement with the expectation from extrapolating 
sideband data into the E\ e f t signal region, 274 ± 19. 

We apply our confidence limit curve to the fitted number of signal events (Fig. EJ) and determine 
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the 90% confidence limit on the signal yield to be 

M f = 22.8, 

which yields a limit on the branching fraction at the 90% confidence level of 

B{B + -> t + v t ) < 4.6 x 10~ 4 at the 90% CL 



0) 
(10) 



where we have assumed equal production of B + B~ and B°B° at the T(4S) resonance. The ex- 
pectation from studies in MC containing zero signal events was given in Sec. El as fif = 15.4; the 
comparison of the expected sensitivity and the measured result is shown in Fig. |H1 
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Figure 8: The nominal sensitivity (15.4), obtained from the median of 10,000 toy MC experiments, 
is compared to the result from data (22.8). 
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Figure 9: The relationship between the confidence level and the branching fraction limit both in 
the case where systematic error is (dashed) and is not (solid) incorporated. Lines are drawn at the 
locations of the 90% confidence limits for both cases. 

After incorporating the systematic uncertainty on the signal efficiency (Fig. |§J) we find an upper 
limit on the branching fraction of 

B{B + ^t + u t ) < 4.9 x 10" 4 at the 90% CL (11) 
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A separate BABAR analysis to search for B + — > t + u t has also been performed where B mesons 
are resonstructed via the semi-exclusive hadronic modes B~ — > D^°X~, where X~ represents a 
combination of up to five charged pions or kaons and up to two ir° candidates. The semi-exclusive 
hadronic reconstruction analysis is statistically independent of the present work and has obtained 
the limit ^S] 

B(B + -» t + u t ) < 7.7 x 1(T 4 at the 90% CL (12) 

The results of this statistically independent analysis are combined with that reported here using 
the same modified frequentist method. For each analysis, 5000 toy MC experiments are generated 
each for a branching fraction hypothesis between 1 x 10~ 4 and 10 x 10 -4 in steps of 2.5 x 10 -6 . The 
number of events in each toy experiment that contains signal and background is denoted n^ s+b y 
For each branching fraction hypothesis a separate sample of 5000 background-only toy experiments 
is also generated; the number of events in a given experiment in this sample is denoted nn,y For a 
given hypothesis i (where i = (s + b), (6), or data) the likelihood estimators, Qi, are calculated in 
each toy experiment: 

* = %^ (13) 
L(ni\b) 

where C{rij\s + b) is the likelihood of a given toy experiment with m events being consistent with 
the (s + b) hypothesis (in equation 01 is set to the number of signal s, corresponding to the 
branching fraction hypothesis, and is set to the number of background b). Systematic errors are 
incorporated by each analysis into the Qi. They are combined as follows: 

gcombined = gSemi-leptonic x ghadronic ^ 

The confidence level for a particular branching ratio hypothesis is given by counting the number 
of experiments, ^(Qi<Q data ), that have a value of Q( s +b) (Q(b)) l ess than that for data (Qdata) and 
taking the ratio 

qj^ (Q(s+b)<Qdata) (15) 



N (Q(b)<Qdata) 



From this a combined BABAR result is obtained: 



B(B + -> T + u T ) combined < 4.1 x 10 -4 at the 90% CL. (16) 



6 Summary 



We have performed a search for the decay process B + — > t + v t . To accomplish this a sample of 
semi-leptonic B decays (B~ — > D°£~i , iX) has been used to reconstruct one of the B mesons and 
the remaining information in the event is searched for evidence of B + — > t + v t . We find no evidence 
for this decay process and set a preliminary limit on its branching fraction of 

B{B + -» t + v t ) < 4.9 x 10~ 4 at the 90% CL 

By combining this analysis with a statistically independent B + — > t + v t search performed using a 
semi-exclusive B reconstruction we find the preliminary combined limit: 

B{B + -» T + u T ) comUn ed < 4.1 x 10~ 4 at the 90% CL 
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